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Abstract 


Dr.  Merrill  K.  King** 
Atlantic  Reserrch  Cor  x>ration 
Alexandria,  Virginia  22314 


Development  of  solid  rocket  motor  designs  which 
Tesult  in  high  velocity  flows  of  product  gases 
across  burning  propellant  surfaces  (notably,  nozzle- 
less rocket  motors)  is  leading  to  increased  occur- 
rence of  erosive  burning.  In  this  paper,  a physi- 
cally realistic  picture  of  the  effect  of  such  cross- 
flows  on  composite  propellant  combustion,  based  on 
the  bending  of  columnar  diffusion  flames  by  the 
crossflow,  is  presented.  This  bending  results  in 
shifting  of  the  diffusion  flame  heat  release  zone 
toward  the  surface,  with  consequent  increased  heat 
feedback  flux  from  this  flame  to  the  surface  and 
thus  increased  burning  rate.  A relatively  simple 
analytical  model  based  on  this  picture  is  developed 
for  prediction  of  propellant  burning  rate  as  a func- 
tion of  pressure  and  crossflow  velocity,  given  only 
zero-crossflow  burning  rate  versus  pressure  data. 
Model  predictions  and  experimental  results  are  com- 
pared, with  reasonably  good  agreement  being  found. 

Introduction 

Requirements  for  ever  higher  propellant  loading 
fractions  in  solid  propellant  rocket  motors  and  for 
higher  thrust-to-weight  ratios  have  led  to  develop- 
ment of  centrally  perforated  grain  configurations 
with  relatively  low  port-to-throat  area  ratios. 

This,  in  turn,  results  in  high  velocities  of  propel- 
lant gases  across  burning  propellant  surfaces  in  the 
aft  portions  of  these  grains,  leading  to  erosive 
burning.  Moreover,  a series  of  studies  has  demon- 
strated that  the  nozzleless  rocket  concept  offers 
significant  economic  advantages  over  a more  conven- 
tional rocket  system  when  considered  for  some  tacti- 
cal weapon  systems.  This  concept  requires  that  the 
flow  within  the  bore  or  central  perforation  of  a 
grain  accelerate  to  the  point  that  sonic  conditions 
are  achieved  at  the  aft  end.  In  this  situation, 
the  high  velocity  environment -leads  to  a realm  of 
erosive  burning  not  previously  considered.  The 
eff~?rs  are  critical  in  that  the  erosive  burn  rate 
contributions  strongly  influence  performance  level, 
performance  repeatability  and  thrust  misalignment. 
More  than  in  any  conventional  motor,  the  exact  ero- 
sive burn  rate  behavior  must  be  held  constant  from 
batch  to  batch  if  reproducibility  is  not  to  be  a 
problem.  The  performance  sensitivity  of  a nozzle- 
less motor  to  erosion  is  due  to  the  fact  that  the 
maximum  erosion  occurs  at  the  choke  point  in  the 
base.  Since  this  point  is  the  effective  throat  area, 
and  the  throat  area  versus  time  is  thus  a function 
of  regression  rate,  the  result  is  a chamber  pressure 
history  which  varies  strongly  with  erosion. 

In  a nozzleless  motor,  two  parameters  which 
affect  burning  rate,  pressure  and  crossflow  velocity, 
vary  strongly  from  the  fore  end  to  the  aft  end  of 
the  grain,  static  pressure  decreasing  and  crossflow 
velocity  increasing  with  distance  from  the  head  end 
of  the  grain.  Assuming  that  an  erosive  burning  rate 
expression  uf  the  form,  r/rQ  * 1 +■  k-jM,  is  appli- 
cable (with  rc  - bpn) , it  may  be  shown*  that  for 
constant  port  area  along  the  grain: 

♦This  effort  was  supported  by  AFOSR  under  Contract 
No.  F44620-76-C-0023 , monitored  by  Maj.  Thomas  Meier 
♦♦Chief  Scientist,  Research  & Technology,  Member  A1AA 


raft/rfore=  (k3+1>/(^  + »"  (1> 

Values  of  (raft/rfore) initial  as  a function  of  the 
erosivity  constant  (IC3)  and  the  burning  rate  expon- 
ent (n)  are  presented  in  Table  I for  y =1.25.  As 
may  be  seen,  for  the  case  of  no  erosion  (IC3  = 0) 
the  aft-end  will  recede  more  slowly  than  the  fore 
end,  due  to  lower  pressure  at  the  aft  end.  As  k3 
increases,  the  rafc'/rfore  ratio  also  increases, 
going  through  unity  (generally  desirable)  at  a value 
of  k3  which  depends  on  the  burning  rate  exponent. 

The  results  of  Table  I give  some  indication  of  the 
sensitivity  of  nozzlelcss  motor  design  to  the  ero- 
sive burning  characteristics  of  the  propellant  and 
thus  point  out  the  importance  of  information  regard- 
ing the  propellant's  erosive  burning  characteris- 
tics to  the  designer. 


Development  of  a better  understanding  of  the 
effects  of  crossflows  on  solid  propellant  combustion 
is  needed  for  accomplishment  of  accurate  motor  per- 
formance predictions  in  terms  of  both  mean  interior 
ballistics  analysis  and  prediction  of  motor  stability 
characteristics.  With  such  understanding,  the  motor 
designer  can  either  design  his  grains  to  compensate 
for  mean  erosive  burning  effects  on  grain  burn  pat- 
tern, or,  knowing  how  propellant  formulation  para- 
meters affect  erosion  sensitivity,  vary  propellant 
parameters  in  such  a way  as  to  optimize  these  effects. 
A review  of  the  literature  indicates  that  there  is 
no  unifying  predictive  model  for  erosive  burning  of 
solid  propellants. 


General  observations  of  importance  from  past  ex- 
perimental studies^~*9ncludc: 

1.  Plots  of  burning  rate  versus  gas  velocity  or 
mass  flux  at  constant  pressure  arc  usually 
not  fitted  best  by  a straight  line. 

2.  Threshold  velocities  and  "negative"  erosion 
rates  are  often  observed. 

3.  Slower  burning  propellants  are  more  strongly 
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affected  by  crossflows  than  higher  burning- 
rate  formulations. 

4.  At  high  pressure,  the  burning  rate  under 
erosive  conditions  tends  to  approach  the 
same  value  for  all  propellants  (at  the  same 
flow  velocity)  regardless  of  the  burning 

I rate  of  the  propellants  at  zero  crossflow. 

5.  Erosive  burning  rates  do  not  depend  upon  gas 
temperature  of  the  crossflow  (determined  from 
tests  in  which  various  "driver  propellants" 
products  are  flowed  across  a given  test  pro- 
pellant) . 

There  is,  however,  very  little  data  available 
for  high  crossflow  velocities  (greater  than  M ~ 0.3). 
In  addition,  there  has  been  no  study  in  which  var- 
ious propellant  parameters  have  been  systematically 
varied  one  at  a time.  Such  a study  is  necessary 
for  determination  of  erosive  burning  mechanisms  and 
proper  modeling  of  the  erosive  burning  phenomena. 
Much  of  the  past  work  has  not  resulted  in  instant- 
aneous (as  opposed  to  averaged  over  a range  of  pres- 
sure and  crossflow  velocity)  measurements  of  eros- 
ive burning  rates  under  well  characterized  local 
flow  conditions. 

From  the  above  disci:--  ion,  it  is  apparent  that 
development  of  an  analvtical  model  of  erosive  burn- 
ing, properly  describing  the  physical  effects  which 
result  in  augmentation  of  solid  composite  propel- 
lant burning  rate  by  crossflows,  coupled  with  an 
experimental  effort  to  systematically  define  the 
effects  of  various  formulation  parameters  on  eros- 
ive burning  at  crossflow  velocities  up  to  Mach  1 is 
of  great  importance  to  the  design  and  development 
of  advanced  solid  rocket  systems. 

Background:  Existing  Models 

The  objectives  of  a theoretical  model  of  erosive 
burning  are  to  provide  a means  of  predicting  the 
sensitivity  of  propellant  combustion  rate  to  gas 
flow  parallel  to  the  ablating  surface  and  to  indi- 
cate what  effect  various  formulation  parameters 
have  on  this  sensitivity.  An  acceptable  model  must 
account  for:  (1)  any  effects  observed  when  cross - 

flow  gas  temperature  is  varied;  (2)  observed  pres- 
sure dependency;  and  (3)  nullification  of  catalyst 
activity  under  erosive  conditions.  This  model 
should  provide  an  explanation  of  the  observed  be- 
havior in  terms  of  the  hydrodynamic  conditions  in- 
duced by  a crossflow  coupled  with  the  chemical  and 
physical  processes  that  constitute  the  propellant 
deflagration  mechanism. 

Over  the  years,  a large  number  of  erosive  burn- 
ing models,  of  varying  degrees  of  sophistication, 
have  been  developed  - a list  of  models  examined  by 
this  author  is  presented  as  Table  II.  These  models 
generally  fall  into  one  of  three  categories,  as  in- 
dicated. The  first  category  is  based  on  the  assump- 
tion that  the  erosive  burning  is  driven  by  increased 
heat  transfer  from  the  mainstream  gas  flow  associat- 
ed with  increased  heat  transfer  coefficient  with 
increased  mass  flux  parallel  to  :he  grain  surface. 

The  best-known  and  most  widely  used  erosive  burning 
model,  that  of  Lenoir  and  Robillard*^  falls  into 
this  category.  In  this  model,  the  authors  calculate 
the  total  burning  rate  (r)  as  the  sum  of  the  normal 
(no  crossflow)  burning  rate  and  a second  erosive 
term  resulting  from  heat  transfer  from  the  "core" 
flow  to  the  propellant  surface.  This  equation  en- 
tails an  a priori  assumption  that  the  pressure-de- 
pendent "base"  rate  (rQ)  is  unaffected  by  an 


increase  in  total  rate  at  a given  pressure,  very  un- 
likely within  the  constraints  of  other  assumptions 
in  the  model.  This  problem  has  been  discussed  in 
detail  by  King^^,  with  derivation  of  a modified 
Lenoir  and  Kobillard  expression  allowing  for  the 
coupling  of  flame  standoff  distance  with  burning 
rate.  While  Lenoir  and  Robillard  assume  r = rn  + re, 
allowance  for  this  coupling  results  in  r = (r0^/r)  + 
re.  In  physical  terms,  Lenoir  and  Robillard  have 
failed  to  account  for  the  fact  that  increased  burn- 
ing rate,  caused  by  erosive  feedback  at  constant 
pressure,  results  in  the  propellant  flame  being 
pushed  further  from  the  surface,  decreasing  its  heat 
feedback  rate,  and  thus  decreasing  the  propellant 
burning. rate  part  of  the  way  back  toward  the  base 
rate. 

Table  II.  General  Types  of  Models  of  Erosive  Burning  Developed  to  Date. 


1.  MODELS  BASED  0.4  HEAT  TRANSFER  FROM  A "CORE  GAS"  IN 
THE  PRESENCE  OF  CROSS  FLOW 

LENOIR  & ROBILLARD  (12) 

BURICK  AND  OSBORN  (13) 

ZUCROW.  OSBORN  AND  MURPHY  (14) 

SADERHOLM  (3) 

MARKLUND  (8) 

JOJIC&  BLAGOJEVIC  (15) 

2.  MOOELS  BASED  ON  ALTERATION  OF  TRANSPOR  T PROPERTIES 
IN  REGION  FROM  SURFACE  TO  FLAME  ZCNF  BY  CHOSSFLOV* 
GENERALLY  DUE  TO  TURBULENCE  EFFECTS  INCLUDES 
EFFECTS  ON  CONDUCTIVITY  FROM  FLAME  ZONE  BACK  TO 
PROPELLANT  ANO  EFFECTS  ON  TIME  FOR  CONSUMPTION  OF 
FUEL  POCKETS  LEAVING  SURFACE. 

SADERHOLM.  BIDDLE.  CA  VENY.  et  al  1 16> 

LENGELLt  (17) 

CORNER  (DOUBLE  BASE)  (18) 

VANDENKERCKHOVE  (DOUBLE  BASE)  (19) 

ZELDOVICH  (DOUBLE  BASE)  (201 
VILYUNOV  (DOUBLE  BASE)  (11) 

GECKLER  (211 

3.  MODELS  BASED  ON  CHEMICALLY  REACTING  BOUNDARY 
LAYER  THEORY  (HOMOGENEOUS  SYSTEMS  ONLY) 

TSUJI  (22) 

BEDDINI.et  *1  (23) 

KUO,  et  •)  124) 


4.  OTHER 


KLIMOV  (25) 
MOLNAR  (26) 
MILLER  (271 
KING 


A more  general  weakness  of  model:  in  the  first 
category  is  that  these  models  predict  substantial 
dependence  on  the  temperature  of  the  core  gas:  such 

dependence  was  found  by  Marklund  and  Lake^  to  be 
completely  absent.  Analysis  of  the  Lenoir  and 
Robillard  treatment  indicates  that  the  erosive  con- 
tribution to  burning  rate  (re)  is  given  by: 


r0.8  0.2  , T v 

G u-  (T  - T ) . 

gas  core  gas  s 


(2) 


for  a given  test  propellant  and  geometry.  But,  at 
fixed  crossflow  velocity  and  pressure,  G is  inverslv 
proportional  to  the  core  gas  (driver  propellant  pro- 
ducts) temperature  while  ugas  is  roughly  directly 
proportional  to  this  temperature.  Therefore, 

r.  « T-0'6  (T  -T  ) . (3) 

e core  gas  core, gas  s 

However,  Marklund  and  Lake  performed  a set  of  exper- 
iments in  which  crossflow  velocity  and  pressure  were 
held  constant  while  the  driver  propellant  was  chanc- 
ed from  a 1700°K  propellant  to  a 2500°K  propellant, 
with  Ts  being  approximately  800°K  in  both  cases. 
Thus,  the  Lenoir  and  Robillard  theory  would  indi- 
cate that: 


-0.6  version  of  the  second-order  closure  method  of 

re  , 2500°K  driver  1700  / 2500  \ = ..  Donaldson  Is  used  to  calculate  the  details  of  the 

r 11Ancv  , . 900  \1700  / = ^ turbulent  flow  field.  To  date,  this  flow  field 

e*  1700  K driver  analysis  has  been  coupled  only  with  a very  simple 

That  is,  wit1'  the  higher  driver  gas  temperature  model  of  propellant  combustion  in  which  the  propcl- 

case,  the  erosive  burning  rate  component  of  the  tot-  lant  surface  temperature  is  assumed  to  be  indcp?nd- 

al  burning  rate  should  be  50  percent  higher  than  ent  of  burning  rate,  the  mass  burning  rate  is  assum- 

that  for  the  low  driver  gas  temperature.  However,  ed  to  be  directly  proportional  to  the  heat  flux 

as  mentioned,  Marklund  and  Lake  observed  no  differ-  from  the  gas  to  the  surface,  and  the  gas  phase  re- 

ence  in  erosive  rates  in  the  two  cases.  This  ob-  action  is  assumed  to  be  described  by  a single  step 

served  lack  of  dependence  of  the  erosive  burning  homogeneous  reaction  which  is  first  order  with  re- 

rate on  core  gas  temperature  tends  to  cast  doubt  on  spect  to  fuel  concentration  and  first  order  with 
all  models  in  the  first  category  of  Table  II.  respect  to  oxidizer  concentration.  Under  these 

conditions,  the  model  predicts  a threshold  velo- 

The  second  category  of  models  listed  in  Table  II  city  for  erosion  followed  by  a quasi-linear  depend- 
includes  models  based  upon  the  alteration  of  trans-  ence  of  burning  rate  on  crossflow  velocity  due  to 

port  properties  in  the  region  between  the  gas  flame  turbulence  entering  the  propellant  flame  zone, 

and  the  propellant  surface  by  the  crossflow,  gener- 
ally due  to  turbulence  effects.  Included  in  this 

category  are  models  in  which  the  thermal  conduct i-  Ktto  and  Razdan^are  also  using  a second -order 

vity  in  this  region  is  raised  by  turbulence  and  turbulence  closure  model  for  characterization  of  the 

models  in  which  the  time  for  consumption  of  fuel  gas  flow  field  in  erosive  burning  situations,  the  clo- 

pockets  leaving  the  surface  is  reduced  by  the  effects  sure  model  being  used  differing  from  that  being 

of  turbulence  on  diffusivity.  Four  of  these  models  used  by  Beddini.  In  addition,  postulated  flame 

were  developed  for  double-base  propellants  as  indi-  mechanisms  (the  details  of  which  are  unknown  to  this 

cated,  and  will  not  be  reviewed  here.  Of  the  com-  author)  for  composite  propellants  are  being  coupled 

posite  propellant  models  in  this  category,  that  of  into  the  analysis.  At  this  time,  the  governing 

Lengelle^7  appears  to  be  the  most  advanced.  The  equations  have  been  developed  and  boundary  condi- 

basic  propellant  combustion  mechanism  assumed  is  the  tions  defined,  but  the  equation  solving  procedure 

granular  diffusion  model  in  which  pockets  of  fuel  has  not  been  completed, 

vapor  leave  the  surface  and  burn  away  in  an  oxidizer 

continuum  at  a rate  strongly  dependent  upon  the  rate  The  Klimov  model  \s  mainly  aimed  at  calculation 

of  micromixing  of  the  oxidizer  vapor  into  the  fuel  of  threshold  crossflow  velocities  (below  which  the 

vapor  pocket.  The  driving  mechanism  by  which  the  propellant  is  unaffected  by  crossflow).  Klimov 

crossflow  is  assumed  to  increase  the  burning  rate  claims  that  the  threshold  velocity  is  the  main 

is  through  increased  turbulence  associated  with  in-  stream  crossflow  velocity  above  which  the  "turbu- 

creasing  crossflow  raising  the  turbulent  diffusivity  lence  front"  subsides  onto  the  propellant  surface, 

in  the  mixing  region  (thus  increasing  the  rate  of  and  presents  boundary  layer  analysis  procedures  for 

mixing  of  the  binder  and  oxidizer  product  gases)  and  calculating  this  threshold  velocity  as  a function 
raising  the  effective  turbulent  thermal  conductivity,  of  the  transpiration  (blowing)  velocity  of  the  gases 
The  increase  in  thermal  conductivity  increases  the  ablating  from  the  propellant  surface.  In  addition, 
heat  transfer  rate  from  the  flame  to  the  surface,  he  postulates  that  negative  erosion  (sometimes  seen 

while  the  increase  in  mixing  rate  just  offsets  the  at  low  crossflow  velocities)  is  due  to  the  "sitrring 

increase  in  gas  velocity  away  from  the  surface,  with  up"  of  cool  streams  of  binder  decomposition  products 
the  result  that  the  flame  offset  distance  remains  over  the  oxidizer  surface,  leading  to  intensif icat- 

constant.  There  are  several  weaknesses  associated  ion  of  their  cooling  effect  and  to  screening  of  heat 

with  the  Lengelle  model:  (l).the  granular  diffusion  feedback  from  the  diffusion  flame, 

flame  model  is  not  physically  realistic;  (2)  the 

ammonium  perchlorate  monopropellant  flame  is  ignored;  Molnar's  model^,  developed  for  homogeneous  pro- 
and  (3)  the  boundary  layer  treatment  used  to  calcu-  pellants  with  a laminar  crossflow,  is  based  on  an 

late  the  dependence  of  the  effective  turbulent  dif-  assumption  (which  does  not  appear  to  this  author  to 

fusivity  and  conductivity  on  the  crossflow  is  un-  be  substantiated)  that  the  lateral  velocity  gradient 

realistic  in  its  use  of  a 1/7 th  power  velocity  law  at  the  burning  surface  governs  erosive  burning, 

all  the  way  from  the  free-stream  to  the  surface.  Miller^7  assumes  that  the  time  for  a unit  of  propel- 

lant to  be  consumed  is  a linear  sum  of  a chemical 


Of  the  three  models  listed  in  the  third  category 
(models  based  on  chemically  reacting  boundary  layer 
theory)  one  is  complete,  while  the  remaining  two  are 
in  development.  The  completed  model,  by  Tsuji^, 
is  unfortunately  not  particularly  useful  due  to  the 
assumption  of  a totally  laminar  boundary  layer  and 
limitation  to  a situation  where  the  free-stream  vel- 
ocity is  proportional  to  the  distance  from  the  head- 
end  of  the  grain.  Other  simplifications  include 
assumption  of  premixed  s toichic ' .etric  fuel  and  oxi- 
dizer (rendering  the  model  inapplicable  to  composite 
propellant  systems)  and  use  of  one-step  global  kin- 
etics. 

21 

In  the  model  of  Beddini,  et  al  , primary  empha- 
sis is  placed  on  analysis  of  a well-developed  tur- 
bulent flow  field  in  a propellant  grain  port  for 
definition  of  turbulent  transport  of  heat,  mass, 
and  momentum  in  the  boundary  layer.  An  extended 


reaction  time  and  "the  time  required  for  turbulent 
transport  of  heat  to  the  propellant  surface"  - such 
an  additivity  approach  does  not  appear  to  this 
author  to  be  physically  realistic. 

Of  the  models  briefly  discussed  above,  those  of 
Lengelle^7,  Beddini,  et  al^3>  and  Kuo  and  Razdan^ 
appear  to  be  the  most  advanced  (although  it  is  not 
clear  at  this  time  how  the  latter  two  teams  will 
build  the  complex  heterogeneous  flame  slruclure 
associated  with  composite  propellant  combustion  into 
their  fluid  dynamic  boundary  layer  treatments) . 

Common  to  all  three  of  these  models  is  the  assump- 
tion that  the  increase  in  propellant  burning  rate 
associated  with  crossflow  results  from  turbulence 
associated  with  this  crossflow  penetrating  between 
the  propellant  gas  flame  zone(s)  and  the  surface, 
causing  increases  in  mass  and  energy  transport  rates. 
However,  for  a typical  propellant  containing  oxidizer 
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particles  with  diameters  of  from  10  to  50  um,  dif- 
fusion flame  offset  distances  may  be  calculated  to 
be  typically  of  the  order  of  one-quarter  to  one- 
half  of  the  particle  diameter  or  2.5  to  25  um.  On 
.the  other  hand,  for  a crossflow  velocity  of  200  m/s 
(650  ft/sec),  the  universal  u+,  y+  correlation 
^(transpiration  effects  neglected)  indicates  a lami- 
nar sublayer  thickness  of  approximately  10  um  and 
a buffer  zone  of  approximately  50  um,  full  turbu- 
lence not  being  achieved  closer  than  60  um  from  the 
propellant  surface.  Moreover,  transpiration  of  the 
binder  and  oxidizer  decomposition  gases  from  the 
propellant  surface  will  tend  to  increase  the  thick- 
ness of  these  zones.  Thus  it  is  not  at  all  certain 
that  crossflow-induced  turbulence  does  penetrate 
into  the  zone  between  the  propellant  surface  and 
the  gas  phase  flame  zone(s).  In  addition,  even  if 
the  turbulent  region  does  extend  into  this  zone, 
in  order  for  the  eddies  to  have  significant  effect 
on  mixing  and  thus  on  heat  and  mass  transfer,  they 
must  be  considerably  smaller  than  the  flame  offset 
distance  - that  is,  they  must  be  on  the  order  of 
1 um  in  diameter  or  less.  It  is  not  at  all  clear 
to  this  author  that  a significant  degree  of  turbu- 
lence of  this  scale  will  be  induced  in  the  zone  be- 
tween the  propellant  surface  and  the  gas  phase  flame 
zone(s)  by  crossflows  even  up  to  Mach  1,  more  than 
an  order  of  magnitude  above  typical  erosive  burning 
threshold  velocities.  Accordingly  an  alternate  pos- 
sible mechanism  for  erosive  burning  of  composite 
propellants  is  presented  below. 


face  in  exactly  the  same  way,  the  variation  will  not 
be  great.  Ihus  fuel  and  oxidizer  gas  columns  leav- 
ing the  surface  will  not  flow  perpendicular  to  the 
surface  (as  they  would  in  the  absence  of  crossflow) 
but  at  an  angle  of  approximately  20  to  25  degrees 
from  parallel  with  the  surface. 

The  important  feature  of  this  picture  is  that 
any  diffusion  flame  at  the  AP-binder  boundaries  is 
bent  over  toward  the  propellant  surface  by  the  cross - 
flow  velocity.  Since  the  deflection  of  this  mixing 
column  or  cone  can  be  shown  to  cause  the  distance 
from  the  base  to  the  tip,  measured* perpendicular  to 
the  base,  to  decrease,  the  height  above  the  propel- 
lant at .which  any  given  fraction  of  the  mixing  of 
AP  products  and  fuel  decomposition  products  is  com- 
plete should,  therefore,  be  decreased  and  the  dis- 
tance from  the  propellant  surface  to  the  "average" 
location  of  the  diffusion  flame  should  also  be  de- 
creased. This,  in  turn,  will  increase  heat  feed- 
back and  thus  increase  burning  rate.  The  schematic 
of  a composite  propellant  erosive  burning  model 
based  upon  this  picture  is  shown  in  Figure  1. 
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Model  Development 

In  development  of  a proper  model  of  erosive  burn- 
ing, it  is  necessary  that  a physical-chemical  mech- 
anism for  the  "normal"  (no  crossflow)  burning  of 
such  propellants  be  specified,  that  the  boundary 
layer  flow  be  properly  described  (theoretically  or 
empirically)  and  that  the  descriptions  of  these  pro- 
cesses be  properly  coupled. 

Considering  first  the  flow  field,  it  is  inform- 
ative to  estimate  flow  profiles  and  angles  near  the 
surface  of  a composite  propellant  for  a typical  er- 
osive burning  situation.  As  an  example,  let  us  ex- 
amine a case  where  the  operating  pressure  is  6.89 
(10^)  N/m^  (1000  psi) , the  propellant  flame  temper- 
ature is  3000  °K,  the  crossflow  mainstream  velocity 
is  200  m/sec  (650  ft/sec),  the  characteristic  length 
dimension  for  determining  Reynold* s Number  is  15  cm 
(0.5  feet)  and  the  propellant  burning  rate  is  1.25 
cm/sec  (0.5  in/sec).  In  this  case,  the  gas  velocity 
away  from  the  surface  r.t  the  flame  temperature  is 
approximately  4 m/sec  (13  ft/sec).  Using  Mickley 
and  Davis^  flow  profile  data  for  boundary  layer 
profiles  in  the  presence  of  transpiration,  we  esti- 
mate that  the  crossflow  velocity  10  um  from  the 
propellant  surface  is  about  10  m/sec  (30  ft/sec). 

A simplified  energy  balance  equating  the  heat  feed- 
back flux  from  a flame  sheet  above  a propellant  sur- 
face to  the  value  required  for  preheating  and  vapor- 
izing the  solid  ingredients  at  a regression  rate  of 
1.25  cm/sec  (0.5  in/sec)  indicates  that  the  gas- 
phase  flame  must  be  on  the  ordei  of  10  um  from  the 
surface.  Thus,  at  the  position  of  the  gas-phase 
flame  front,  the  velocity  component  away  from  the 
propellant  is  about  4 m/sec  while  the  velocity  com- 
ponent parallel  to  the  surface  is  10  m/sec,  and  the 
resultant  flow  vector  makes  an  angle  with  the  pro- 
pellant surface  of  only  22  degrees.  While  this  vec- 
tor will  vary  with  distance  from  the  surface,  since 
the  velocity  components  normal  to  and  parallel  to 
the  surface  do  not  scale  with  distance  from  the  sur- 
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Figure  1.  Schematic  of  Geometrical  Model  of  Erosive  Burning 
(Two-Flame  Model). 

In  the  first  part  of  the  figure,  the  flame  pro- 
cesses occurring  in  the  absence  of  crossflow  are 
depicted.  There  are  two  flames  considered,  an  ammon- 
ium perchlorate  deflagration  monopropellant  flame 
close  to  the  surface  and  a columnar  diffusion  flame 
resulting  from  mixing  and  combustior  of  the  AP  de- 
flagration products  and  fuel  binder  pyrolysis  pro- 
ducts at  an  average  distance  somewhat  further  from 
the  surface.  Three  important  distance  parameters 
considered  are  the  distance  from  the  propellant  sur- 
face to  the  "average"  location  of  the  kinetically 
controlled  AP  monopropellant  heat  release  (Lj) , the 
distance  associated  with  mixing  of  the  oxidizer  and 
fuel  for  the  diffusion  flame  (L^f  f)  , and  the  dis- 
tance associated  with  the  fuel-oxidizer  reaction 
time  subsequent  to  mixing  (I-Kin)  • A heat  balance 
between  heat  feedback  from  these  two  flames  and  the 
energy  requirements  for  heating  the  propellant  from 
its  initial  temperature  to  the  burning  surface  tem- 
perature and  decomposing  it  yields  (assuming  that 
the  heat  feedback  required  per  unit  mass  of  propel- 
lant consumed  is  independent  of  burning  rate) : 
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The  situation  pictured  as  prevailing  with  a cross 
flow  is  shown  in  the  second  part  of  Figure  1.  Since 
Lj  and  are  both  kinetically  controlled  and  are 

thus  simply  proportional  to  a characteristic  react- 
ion time  (which  is  assumed  to  be  unaffected  by  the 
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crossflow)  multiplied  by  the  propellant  gas  velocity 
normal  to  the  surface  (which  for  a given  formulation 
is  fixed  by  burning  rate  and  pressure  alone)  these 
distances  are  fixed  for  a given  formulation  at  a 
given  burning  rate  and  pressure,  independent  of  the 
crossflow  velocity.  Of  course,  since  crossflow  vel- 
ocity affects  burning  rate  at  a given  pressure  thro- 
ugh its  influence  on  the  diffusion  process  as  dis- 
cussed below,  Lj  and  I^n  are  influenced  through  the 
change  in  burning  rate,  but  this  is  simply  coupled 
into  a model  by  expressing  Lj  and  Lj^n  as  explicit 
functions  of  burning  rate  and  pressure  in  that  model. 
The  important  point  is  that  they  can  be  expressed 
as  functions  of  these  two  parameters  alone  for  a 
given  propellant.  However,  the  distance  of  the  mix- 
ing zone  from  the  propellant  surface  is  directly 
affected  by  the  crossflow.  It  may  be  shown  through 
geometrical  arguments  coupled  with  the  columnar  dif- 
fusion flame  height  analysis  presented  by  Schultz, 
Penner  and  Green^)  # that  Lpiff  measured  along  a 
vector  coincident  with  the  resultant  crossflow  and 
transpiration  velocities  should  be  approximately  the 
same  as  Lpiff  in  the  absence  of  a crossflow  at  the 
same  burning  rate  and  pressure  (except  at  very  high 
ratios  of  local  crossflow  velocity  to  transpiration 
velocity).  That  is,  the  magnitude  of  Lp^ff  is  essen- 
tially independent  of  the  crossflow  velocity,  altho- 
ugh its  orientation  is  not.  Thus,  the  distance  from 
the  surface  to  the  "average"  mixed  region  is  decre- 
ased to  sin  0 where  0 represents  the  angle 

between  the  surface  and  the  average  flow  vector  in 
the  mixing  region.  The  heat  balance  at  the  propel- 
lant surface  now  yields: 
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for  burning  in  the  absence  of  crossflow.  A 
regression  analysis  using  no-crossflow  burn- 
ing rate  data  is  performed  to  obtain  best  fit 
values  for  A3,  A4,  and  A5.  (dp  is  the  aver- 
age ammonium  perchlorate  particle  size.  For 
a given  propellant,  the  burning  rate  data 
may  be  just  as  effectively  regressed  on  A3, 
A4,  and  A5<iS,  eliminating  the  necessity  of 
actually  defining  an  effective  average  parti- 
cle size.) 


From  these  results,  expressions  are  obtained 
for  Lj,  Lpiff,  and  Lgin  as  functions  of  burn- 
ing rate  (or  m)  and  pressure. 

These  expressions  are  combined  with  an  analy- 
sis of  the  boundary  layer  flow  (which  gives 
the  crossflow  velocity  as  a function  of  dis- 
tance from  the  propellant  surface,  mainstream 
velocity,  and  propellant  burning  rate)  to 
permit  calculation  of  the  angle  9 (Figure  1) , 
LI>  LDiff»  Tfcin*  ant*  f°r  a given  pressure 
and  crossflow  velocity. 


In  the  derivation  of  a burning  rate  expression 
for  a composite  propellant  in  the  absence  of  a cross- 
flow,  an  energy  balance  at  the  propellant  surface 
is  first  written  as:  (See  Figure  1.) 
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This  picture  has  been  used  as  the  basis  of  devel- 
opment of  a model  for  prediction  of  burning-rate 
versus  pressure  curves  at  various  crossflow  velo- 
cities, given  only  a curve  of  burning  rate  versus 
pressure  in  the  absence  of  crossflow.  This  model 
employs  no  empirical  constants  other  than  those 
obtained  from  regression  analysis  of  the  no-cross- 
flow  burning  rate  data.  Thus,  although  it  is  not 
as  powerful  as  a model  which  would  permit  prediction 
of  erosive  burning  phenomena  with  no  burning  rate 
data  at  all,  but  only  propellant  composition  and 
ingredient  size  data,  it  is  still  a very  useful  tool 
in  that  it  permits  prediction  of  erosive  burning 
characteristics  given  only  relatively  easily  obtain- 
ed strand-bomb  burning  rate  data.  (By  comparison, 
the  Lenoir  and  Robillard  model  employs  two  free  con- 
stants which  are  adjusted  to  provide  a best  fit  of 
erosive  burning  data  for  a given  propellant  and 
since  since  these  constants  vary  from  propellant  to 
propellant,  the  Lenoir  and  Robillard  model  does  not 
permit  a priori  erosive  burning  predictions  for  new 
propellants  without  some  erosive  burning  data,  where 
as  the  model  presented  here  does  not  require  such 
data.) 

The  general  approach  followed  in  development  of 
this  model  is: 

1.  The  expressions  for  Lj»  and  LKin  as 

functions  of  burning  rate  (or  burning  mass 
flux,  m) , pressure,  and  propellant  properties 
are  derived  and  substituted  into  a propellant 
surface  heat  balance. 

2.  The  resulting  equation  is  worked  into  the 
form  (developed  in  succeeding  paragraphs): 


The  first  term  of  this  equation  represents  heat  flux 
from  the  final  flame  to  the  surface,  the  second  re- 
presents heat  flux  from  the  AP  monopropellant  flame, 
and  the  third  represents  the  heat  flux  requirements 
for  ablation  of  the  propellant  at  the  mass  flux,  m. 
Several  simplifying  assumptions  are  obviously  invol- 
ved in  writing  of  the  equation  in  this  form.  Pro- 
bably the  most  important  and  tenuous  of  these  is 
the  assumption  that  Qrx  is  independent  of  burning 
rate  (or  A)  and  of  pressure.  In  the  Zeldovich  pic- 
ture of  solid  propellant  combustion,  where  subsur- 
face exothermic  reactions  with  fairly  high  activa- 
tion energies  are  considered  to  dominate,  this  would 
be  a very  poor  assumption,  but  in  the  generally 
accepted  picture  of  solid  propellant  combustion  in 
this  country,  it  is  not  a bad  approximation.  In 
addition,  it  is  assumed  that  the  surface  temperature 
is  nearly  constant  with  respect  to  pressure  and  burn- 
ing rate,  with  the  resultant  uncoupling  of  this  heat 
balance  equation  from  a surface  regression  rate 
Arrenhius  expression.  Finally,  it  is  assumed  that 
for  the  diffusion  flame,  a distance  associated  with 
mixing  may  be  added  linearly  to  a distance  associat- 
ed with  reaction  delay  to  yield  a total  flame  offset 
distance,  a fairly  gross  simplification. 

The  monopropellant  AP  flame  offset  distance,  Lj , 
may  be  expressed  as  the  product  of  a characteristic 
reaction  time,  Tj,  and  the  linear  velocity  of  gases 
leaving  the  propellant  surface: 
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For  a second-order  gas-phase  reaction  (generally 
assumed),  Tj  is  inversely  proportional  to  pressure, 
and  for  a given  formulation,  the  gas  density  is 


_ 


- 


F or  8';  » 11  sheet 

T jf,*  page  r.  enroll  — 6'i 


directly  proportional  to  pressure,  yielding: 

4 = KjA/P2  (10) 

A similar  analysis  for  l^n  yields: 

4i„  = K2m/P2  (ID 

For  a columnar  diffusion  flame,  it  may  easily  be 
shown^  that  the  diffusion  cone  height,  Lj^f,  may 
be  expressed  as: 

huff  “ (12) 

Equations  8 and  10  - 12  may  be  combined  to  yield: 
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Burning  rate  versus  pressure  data  for  a given  pro- 
pellant in  the  absence  of  a crossflow  may  then  be 
analyzed  via  a fairly  complicated  regression  analy- 
sis procedure  to  yield  values  of  the  constants  A^, 
A^,  and  A5  (or  A5d£)  for  that  given  propellant. 

The  constants  Kj_,  K2,  and  K3  are  related  to  these 
constants  in  turn  by: 
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In  this  analysis  a rough  estimate  of  A2  has  been 
made  to  permit  calculation  of  values  of  K^,  and 
K3  from  the  best-fit  values  of  A3,  A4,  and  A5.  It 
should  be  pointed  out,  however,  that  the  subsequent 
calculations  of  burning  rates. in  crossflows  are  not 
strongly  affected  by  the  estimate  of  A2,  since  the 
same  value  of  A2  is  used  in  that  analysis,  and  thus 
its  effects  essentially  cancel.  The  value  used  for 
most  cases  (except  those  cases  run  to  test  the 
effect  of  A2)  was  2*10”  gm  sec  °K/cm^, 

Data  of  Mickley  and  Davis^  were  used  to  develop 
empirical  expressions  for  the  local  crossflow  vel- 
ocity as  a function  of  distance  from  the  propellant 
surface,  mainstream  crossflow  velocity,  and  trans- 
piration rate  (gas  velocity  normal  to  the  propellant 
surface).  In  this  analysis,  it  was  decided  that  the 
transpiration  velocity  should  be  calculated  as  the 
gas  velocity  normal  to  the  surface  at  the  final 
flame  temperature,  (Mickley  and  Davis  correlations 
are  based  upon  the  ratio  of  mainstream  velocity  to 
transpiration  velocity.)  The  procedure  used  is  out- 
lined in  Table  III. 

The  above  analyses  were  used  in  the  derivation  of 
the  following  eight  equations  in  eight  unknowns  for 
the  burning  of  a given  composite  propellant  at  a 
given  pressure  and  crossflow  velocity: 


= U*  (See  Table  III) 

and  the  function  f of  Equation  24  is  given  in 
Table  III. 

Table  III.  Calculation  of  Cross-flow  Velocity  Profile 
in  Current  Erosive  Burning  Model. 


1.  NEGLECTING  TRANSPIRATION  EFFECTS 
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Implicit  in  Equation  25  is  the  assumption  that 
the  transpiration  velocity  and  the  crossflow  velo- 
city maintain  a constant  ratio  from  very  near  the 
surface  out  to  the  end  of  the  diffusion  zone:  that 

is,  that  the  vector  resultant  is  a straight  line. 
This  approximation  is  probably  not  seriously  in 
error,  and  should  not  strongly  affect  the  results 
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of  the  calculations. 

As  may  be  seen,  the  quantity  A2  appears  in  the 
denominator  of  k|,  K7,  Kj,  K5,  and  Kg.  Thus,  as 
indicated  earlier,  the  effect  of  A2  in  Equation  18 
cancels  out  and  the  predicted  burning  rato  is  de- 
pendent upon  this  parameter  only  to  the  extent  that 
it  affects  the  calculation  of  the  crossflow  velocity 
at  distance  sin  9 from  the  surface.  Paramet- 

ric calculations  with  various  values  of  A2  indicate 
that  this  effect  is  very  weak.  A computer  code  has 
been  developed  to  solve  these  equations  simultan- 
eously, yielding  a predicted  burning  rate  for  a 
given  pressure,  crossflow  velocity  and  set  of  con- 
stants A3,  Ag,  and  Asdp  obtained  from  regression 
analysis  of  no-crossflow  data. 

Comparison  of  Predictions  With  Data 

Original  testing  of  the  model  was  carried  out 
using  a systematic  erosive  burning  data  set  taken 
by  Saderholm-*.  (This  was  the  only  systematic  data 
set  found  in  the  literature  with  sufficient  zero 
crossflow  data  to  permit  evaluation  of  A3,  Ag,  and 
Agd^.)  The  computer  code  described  above  was  used 
to  calculate  burning  rate  versus  pressure  curves  for 
several  crossflow  velocities  studied  by  Saderholm, 
with  and  without  correction  of  the  boundary  layer 
profiles  for  transpiration  effects.  The  results 
are  shown  in  Figures  2 and  3.  As  shown  in  Figure  2, 
neglecting  correction  for  the  effects  of  transpira- 
tion on  the  boundary  layer  profile  results  in  ser- 
ious overprediction  of  the  burning  rates.  However, 
as  shown  in  Figure  3,  agreement  between  predictions 
and  data  is  excellent  when  the  transpiration  cor- 
rection factor  is  included. 

In  parallel  with  this  modeling  effort,  Atlantic 
Research  is  carrying  out  an  experimental  test  pro- 
gram to  obtain  systematic  erosive  burning  data  for 
a series  of  propellant  formulations.  This  experi- 
mental program  is  described  in  some  detail  in  Ref. 

1 and  31.  At  this  time,  a fairly  complete  set  of 
data  covering  a pressure  range  of  10**  to  5 (10°) 
n/m^  (10  to  50  atmospheres)  and  a crossflow  velo- 
city range  of  200  to  700  m/sec  (600  to  2200  ft/sec) 
have  been  obtained  for  one  formulation,  designated 
Formulation  4525.  This  is  a "scholastic"  formula- 
tion, containing  unimodal  ammonium  perchlorate.  The 
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Figure  3.  Erosive  Burning  Mode)  Predictions  and  Comparisons  with 
Saderholm  Data.  Transpiration  Lffects  Included. 

formulation  consists  of  73  weight  percent  20  micron 
diameter  ammonium  perchlorate  and  27  weight  percent 
hydroxy  terminated  polybutadiene  (HTPB)  binder,  with 
a trace  of  carbon  black  added  to  opacify  the  propel- 
lant. Experimental  and  theoretical  results  are  pre- 
sented in  Figures  4 and  5.  As  may  be  seen,  agree- 
ment between  prediction  and  data,  while  not  as  good 
as  with  the  Saderholm  propellant,  is  not  bad.  The 
predicted  curves  for  burning  rate  versus  pressure 
at  various  crossflow  velocities  (Figure  4)  do  seem 
to  group  more  tightly  than  the  data.  That  is,  as 
shown  more  clearly  in  Figure  5,  the  model  tends  to 
overpredict  the  burning  rate  at  low  crossflow  vel- 
ocities and  underpredict  it  at  high  velocities. 
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Figure  4.  Burning-Rate  Versus  Pressure  Data  and  Predictions  lor  Various 
Crossflow  Velocities  for  Formulation  4525  (73/27  AP/HTPB. 
2Q/im  AP). 
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Figure  2.  Erosive  Burning  Model  Predictions  and  Comparisons  with 
Saderholm  Data.  Transpiration  Effects  not  Included. 
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Figure  5.  Burning-Rate  Versus  Crossflow  Data  and  Predictions  for  Various 
Pressures  for  Formulation  4525  (73/27  AP/HTPB,  20/im  AP). 

Summary 

Past  modeling  efforts  in  the  area  of  erosive 
burning  of  solid  propellants  have  been  reviewed  and 
lack  of  a model  which  incorporates  a realistic  des- 
cription of  composite  propellant  combustion  has  been 
noted.  A possible  physical  mechanism  by  which  cross 
flows  may  affect  the  combustion  of  a composite  pro- 
pellant has  been  postulated  and  a mathematical  model 
for  prediction  of  the  burning  rate  of  a composite 
propellant  in  such  a crossflow,  given  only  the  no- 
crossflow burning  rate  versus  pressure  characteris- 
tics of  the  propellant,  has  been  developed.  This 
model  has  been  used  to  predict  remarkably  well  the 
erosive  burning  characteristics  of  a propellant 
studied  by  Saderholm.  In  addition,  reasonable 
agreement  between  predictions  and  data  has  been 
obtained  for  a formulation  recently  characterized 
in  our  test  facility.  Additional  propellants  are 
currently  being  studied,  and  erosive  burning  rate 
data  and  predictions  for  these  formulations  will  be 
compared  for  these  formulations  for  further  testing 
of  the  model  in  the  near  future.  In  addition,  a 
second  generation  model  based  on  this  same  picture 
which  does  not  require  no-crossflow  burning  rate 
versus  pressure  data  but  instead  uses  only  propel- 
lant composition  and  ingredient  size  as  input  is 
currently  under  development. 
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A.  constant  given  by  Equation  17 
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empirical  constants  relating  zero  cross- 
flow burning  rate  to  pressure,  obtained 
by  regression  analysis  of  data 

pre-exponential  in  Vielle's  burning  rate 
law 

average  propellant  heat  capacity 
oxidizer  particle  diameter 

characteristic  length  for  calculation 
of  Reynold's  Number 

friction  factor 

crossflow  mass  flux 

constants  defined  by  Equation  5 

erosivity  constant  in  expression  of 
form,  r/r0=l+K3M 


core  gas 
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^Crossflow 
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constants  relating  standoff  distances 
to  pressure  and  burning  mass  flux, 
given  by  Equations  14-16 

grouped  constants  defined  by  Eqns. 

26  - 33 

diffusion  distance  (Figure  1) 

distance  associated  with  oxidizer/ 
fuel  reaction  subsequent  to  mixing 
(Figure  1) 

distance  associated  with  oxidizer 
monopropellant  reaction  (Figure  1) 

crossflow  Mach  Number 

molecular  weight  of  propellant  pro- 
duct gases 

propellant  burning  mass  flux  (linear 
burning  rate  x propellant  density) 

exponent  in  Vielle's  burning  rate  law 
Pressure 

heat  feedback  flux  from  gas  flames 
to  propellant  surface 

heat  per  unit  mass  involved  in  var- 
ious endothermic  processes  at  or 
below  the  propellant  surface,  e.g., 
binder  pyrolysis  or  AP  sublimation. 

heat  per  unit  mass  involved  in  var- 
ious exothermic  processes  at  or  below 
the  propellant  surface. 

gas  law  constant 
linear  burning  rate 

linear  burning  rate  at  zero  crossflow 

erosive  contribution  to  linear  burn- 
ing rate  (r  - r ) 
o 

linear  burning  rate  of  aft  end  of 
a grain  port 

linear  burning  rate  at  fore  end  of 
a grain  port 

Reynolds  Number 

ammonium  perchlorate  flame  tempera- 
ture 

propellant  flame  temperature 
propellant  surface  temperature 
propellant  bulk  temperature 

temperature  of  core  gas  flowing  past 
a propellant  surface 

friction  velocity  (shear  velocity) 
of  crossflow 

mainstream  crossflow  velocity 

local  crossflow  velocity  at  a given 
distance  from  the  propellant  surface 

blowing  velocity  of  gases  produced 
by  propellant  combustion  normal  to 
the  surface,  evaluated  at  the  final 
flame  temperature 

distance  from  propellant  surface 

dimensionless  distance,  o U*Y/u 

’ “gas  gas 

specific  heat  ratio  of  propellant 
product  gases 
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gas  thermal  conductivity,  with  an  area 
ratio  term  for  each  flame  included 

propellant  product  gas  density 

solid  propellant  density 

angle  between  propellant  surface  and  vector 
resultant  of  crossflow  and  transpiration 
velocities.  (Figure  1) 

characteristic  reaction  time  for  AP  mono- 
propellant reaction 

propellant  product  gas  viscosity. 
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